Cryptoferromagnetic state in superconductor-ferromagnet multilayers 
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We study a possibility of a non-homogeneous magnetic 
order (cryptoferromagnetic state) in heterostructures consist- 
ing of a bulk superconductor and a ferromagnetic thin layer 
that can be due to the influence of the superconductor. The 
exchange field in the ferromagnet may be strong and exceed 
the inverse mean free time. A new approach based on solv- 
ing the Eilenberger equations in the ferromagnet and the Us- 
adel equations in the superconductor is developed. We derive 
a phase diagram between the cryptoferromagnetic and ferro- 
magnetic states and discuss the possibility of an experimental 
observation of the CF state in different materials. 
PACS: 74.80.Dm,74.50.+r, 75.10.-b 

In the last years, the interest in experiments on 
superconducting-ferromagnet {S/F) hybrid structures 
has grown rapidly. Such structures show the coexistence 
of these two antagonistic orderings but their mutual in- 
fluence is still a controversial point In these ex- 
periments, the multilayers contained strong ferromagnets 
like Fe or Gd with the Curie temperature up to 1000 A" 
and superconductors with transition temperatures not 
exceeding 10AT, like Nb or V. 

Naturally, in most theoretical works only the influence 
of the ferromagnet on the superconductivity of S/F sys- 
tems was considered One may argue that a modi- 
fication of the magnetic ordering would need energies of 
the order of the Curie, which is much larger than the su- 
perconducting transition temperature T c . Therefore, any 
change of the ferromagnetic order would be less energet- 
ically favorable than the destruction of the superconduc- 
tivity in the vicinity of the ferromagnet. 

This simple argument was questioned in a recent ex- 
perimental work [0, where Nb/ Fe bilayers were studied 
using different experimental techniques. Direct measure- 
ments using the ferromagnetic resonance showed that in 
several samples with thin ferromagnetic layers (10— 15A) 
the average magnetic moment started to decay at the 
superconducting transition temperature T c . The mea- 
surements were possible only in a limited range of the 
temperatures below T c and the decrease of the magnetic 
moment in this interval reached 10% without any sign 
of a saturation. As a possible explanation of the effect, 
it was assumed in Ref. [ffo| that the superconductivity 
affected the magnetic order causing a domain-like struc- 
ture. 

A possibility of a domain-like magnetic structure in 
presence of superconductivity has been first suggested 



by Anderson and Suhl long ago |Tl|. They argued that 
a weak ferromagnetism of localized electrons should not 
destroy the superconductivity in the conduction band. 
Instead, it may become more favorable energetically to 
build a domain structure called cryptoferromagnetic state 
]i"l| . Later this state was investigated both theoretically 
and experimentally in detail (for review see, e.g. flT2|)- 

In this paper, we investigate theoretically the possi- 
bility of a cryptoferromagnetic-like (CF) state in S/F 
bilayers with parameters corresponding to the structures 
used in the experiments p]-|q,^0|. Such a study is very 
important because it may allow to clarify the question 
about the cryptoferromagnetic state in the experiment 
JTo| and to make predictions for other S/F multilayers. 
From the theoretical point of view, large magnetic en- 
ergies involved make the problem quite non-trivial and 
demand development of new approaches. 

To the best of our knowledge the possibility of a non- 
homogeneous magnetic order in multilayers was consid- 
ered only in Ref. (l3). However, although the authors 
of Ref. jl3) came to the conclusion that the domain-like 
structure due to the interaction with the superconduc- 
tor was possible, the results obtained can hardly be used 
for quantitative estimates. For example, they assumed 
that the period of the structure b had to be not only 
much smaller than the size of the Cooper pair £, but also 
than £y/T c /h, where h is the energy of interaction of 
conduction electrons (CEs) with the localized magnetic 
moments (LMs). In addition, a very rough boundary 
condition at the S/F boundary was used. 

In contrast, we present here a microscopic derivation 
of the phase diagram valid for realistic parameters of the 
problem involved. We will show that the phase transition 
between the CF and ferromagnetic (F) phases is contin- 
uous and the period of the structure b goes to infinity at 
the critical point. The only restrictions we use are 



d <C £f = vo/h, T c <C h -C e 



(1) 



where d is the thickness of the ferromagnetic layer, vq 
and £o are the Fermi-velocity and Fermi-energy. 

Even in the such strong ferromagnet as iron, £p is of 
the order 10A. For weaker ferromagnets like Gd, £f is 
considerably larger and the inequalities ([!]) can be ful- 
filled rather easily. 

We assume that the superconductor occupies the half- 
space x > while the ferromagnetic film occupies the 
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region — d < x < and write the Hamiltonian as 



H=H BCS +7/ dr*+(r) [h(r)a] Q/3 ^(r) + ff 



(2) 



where Hbcs is the usual BCS Hamiltonian (in the pres- 
ence of non-magnetic impurities) describing the super- 
conducting state in the S layer, 7 is a constant which 
will be put to 1 at the end. The second term in Eq. (||) 
stands for the interaction between the LMs of the ferro- 
magnet and the CEs, where h is the exchange field and 
a is the vector containing the Pauli matrices as compo- 
nents. We neglect the influence of the LMs on the orbital 
motion of the CEs since the exchange interaction is the 
dominant Cooper pair breaking mechanism |12| for the 
problem involved. The term Hm describes the interac- 
tion between the LM in the ferromagnet. 

Our aim is to obtain an expression for the free energy 
of the system for different magnetic structures in the F 
layer. To determine the contribution of an inhomoge- 
neous alignment of magnetic spins to the total energy 
we use the limit of a continuous material and replace the 
spins by classical vectors. We assume that the anisotropy 
energy of the ferromagnet is smaller than the exchange 
energy and hence there is no easy axis of magnetization. 
This can definitely be a good approximation for iron with 
a cubic lattice used in the work |Io| . The energy Hm of 
a non-homogeneous structure can be written in the con- 
tinuum limit as 



H 



M 



J {VS x f + (VS v f + {VS z f dV, (3) 



where the magnetic stiffness J characterizes the strength 
of the coupling between LMs in the F layer and Sj's 
are the components of a unit vector. Writing S = 
(0, — sin O, cos 0) and minimizing the energy Hm we ob- 
tain the equation A0 — 0. We consider only the solutions 
of this equation that are of interest for us: 



practice, this is difficult and we simplify the problem con- 
sidering the "dirty limit" I <ti !;oi where I is the mean free 
path and £0 = v/T c is the coherence length of the super- 
conductor in the clean limit, which allows to use the more 
simple Usadel equations |l5). If we assume that |t| <C 1, 
t = (T — T c ) /T c , the Usadel equations together with 
the self-consistency equation can be further reduced to 
the Ginzburg-Landau (GL) equation [[l6|-|l§|]. However, 
the latter equation can be used only sufficiently far from 
the S/F boundary at distances exceeding £ ~ y/^ol. At 
the distances of the order of £ one should write again 
the Usadel equations but in the limit |r| <C 1 they can 
be linearized. This is a conventional scheme of calcula- 
tion for interfaces between superconductors and normal 
metals or ferromagncts. 

Writing the Usadel equations in the ferromagnet may 
not be a good approximation because the exchange en- 
ergy h in realistic cases is not necessarily smaller than 
1/rtr, where T% r the mean free time, and so one should 
write in this region the Eilenberger equations. At the 
end one should match the solutions of all the equations. 

Now we start the calculations following this program. 
The loss of the superconducting energy due to the sup- 
pression of the superconductivity in the S'-layer can be 
found from the solution of the GL equation for the order 
parameter A(r). At distances the proper solution 

is MM 



A(x) = A(T)tanh 



(V2£(T) 



C 



(6) 



where A(T) = 



87T 2 

7C(3) 



t\T c = Aqt 1 / 2 is the value of 



a) e = 0, 



b) 6 = Qy 



(4) 



the order parameter in the bulk superconductor, £(T) = 
\J~^-\ T \~ 1 ^ 2 is the characteristic scale of the spacial vari- 
ation of A (r) , D is the diffusion coefficient in the super- 
conductor, and C is a constant. Substituting A (a:), Eq. 
(^), into the GL free energy functional one can evalu- 
ate the loss of the superconducting energy at the F/S 
interface per unit surface area as function of C fll7| 



The solution a) in Eq. (J|) corresponds to the F state, 
whereas the solution b) describes a CF state with a ho- 
mogeneously rotating magnetic moment. The wave vec- 
tor of this rotation is denoted by Q. The magnetization 
is chosen to be parallel to the FS interface, i.e. to the 
yz— plane. This allows to neglect Meissner currents in the 
superconductor. With all this assumptions the magnetic 
energy Qm (per unit surface area) is given by 



9. 



JdQ 2 



(5) 



The corresponding energy of the F state equals zero. 

The superconducting part of the energy can be calcu- 
lated deriving from Eq. (|^) proper Eilenberger equations 
[pj for the superconductor and the ferromagnet, solv- 
ing these equations and then matching the solutions. In 



6V2 



|t| 3 / 2 (2 + K) (1 - Kf 



(7) 



where K = tanh C . The influence of the ferromagnet on 
the superconductivity is determined by the parameter K 
that will be found by minimizing the total energy. 

The contribution SIm/s of the second term in (||) to the 
total energy has still to be determined. First, we write 
the Eilenberger equation for the magnetic moment h (r) 
depending on coordinates. Introducing the quasiclassical 
matrix Green function g w (r,po) 
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one derives in the standard way the Eilenberger equation 
in the spin®particle-hole space 



2 



[{ wf 3 - i A + ijV + it unp },§] + v V r .g = . (8) 

where po and vo are the momentum and velocity at the 
Fermi-surface. 

In Eq. (|J), fj, i = 1,2,3, are Pauli matrices in the 
particle-hole space, A=fi <S)ia y A(r), V = Re(h(r)er) ® 
1 + Im(h(r)cr) ® T3, and A should be determined self- 
consistently 



A(r) = -~twA T^ < /i2(r,Po,^„) > . 



(9) 



where T] = Ujf /vq and /o is the zero harmonics of the 
function / in the superconductor. When deriving Eq. 
( |l4| ) we used the fact that the Usadel equation is ap- 
plicable in the S- layer at distances down to the mean 
free path I and extrapolated its solution to the interface. 
Only first two spherical harmonics /W f=a /o + vofi were 
kept in the derivation. 

The linearized Usadel equation for the superconductor 
can be written in the standard form 



DV 2 f Q -2\io\f -2A{x)a y = 



where (...)o denotes averaging over the Fermi velocity 
and Ao is the constant of the electron-electron interac- 
tion, v is the density of states. We assume for simplicity 
that Ao = and hence A = in the ferromagnet. At the 
same time, h = in the superconductor. The term i£ imp 
describes scattering by impurities. For a short range in- 
teraction, £™ D = 



(15) 



The general solution of Eq. (|15|) with the boundary 
condition, Eq. (|l4|), and V from Eq. ([l|) can be written 
as 



/o(r,w) = a u (x)a x e 



(16) 



£ (g) . Eq. (g) is complemented by where ^ = ^ exp 
ndition a = 1. Once we know a. we 



the normalization condition g 2 = 1. Once we know g, we 
can determine £Im/s using the expression ]TF 



exp 



Eq. ([!(]) is applicable at 



M/S 



(10) 



Near T c , the anomalous functions / and / + are small and 
g ~ sgn (u>). Then, in the limit T c <gC h the off-diagonal 
component (1,2) in particle-hole space of the equation (^) 
in the region — d < x < is 



vfW + i$Wv*- £5!iM(/CF)_ < /(f) >} 



v V/=-i 



V" = h(x)a z exp(iQya x ) 

h is the strength of the exchange field in the F-layer. 

Assuming that d ^ v /h we can relate the values of 
the function /( F )(vo,r) at the interface, i.e. at x = 0~ 
to the values at the boundary to the vacuum at x = — d 
using the Taylor expansion: 



distances much smaller than £ (T), where the solution 
for A can be approximated by a linear function. One 
can check using the self-consistency Eq. (^) that the 
relative correction to A coming from the exponentially 

decaying part of Eq. (|l6|), is of the order ^ln ^p.^ , 

where lod is the Debye frequency, and we neglected it. 
The coefficients C u and B u can be now determined from 
Eq. (pT[). Using the condition g 2 — 1 and Eq. ( |l(i| ) we 
can find the energy Introducing the dimensionless 

parameters: 



a 2 = 



2h 2 d 2 




v ,r - r d ) ps / (F) (v ,r ) - dd x f {F) {-v ,r ) , 



(12) 



where r = (0,y,z) and r d = (-d,y,z) . Applying gen- 
eral boundary conditions ]l9t to the problem involved we 
conclude that for a perfectly transparent interface the 
function / is continuous at the interface. At the bound- 
ary with the vacuum (x = —d) the function / satisfies 



DT c g 2 ' 2T, 
and using Eq. (^|) one obtains 

&M/S = ^F 3/2j1 K 2 \t\ + V2F 2A K (1 - K 2 ) \t\ 
+*- 1 F 5 /2,i (l-K 2 ) 2 \r\ 2 , 

F m> t = i] 3/2 _ m 2J q„ m V a «K + 9 2 ) + ( 

n>0 



3/2 



(18) 



/^(^,r -r d )=/^(- Ux) r -r d ) 



(13) 



Using Eqs. (|Tl, 12, 1^) and the continuity of / at r = r 
the problem is reduced to the solving of the Usadel equa- 
tion in the superconductor with the following effective 
boundary condition at the interface 



r ] D{d x + d^)/ (r ) + isgn (u) [~Vf + foV* |=0 (14) 



where a n = ir(2n + 1) and vp is the density of states in 
the ferromagnet. The total energy is given by 57 = J7m + 
Qs + fW/Si Eqs. (||, 0, |l^) and is a functions of two 
parameters, K and q, that should be determined from 
the conditions dfl/dK = dfl/dq = 0. The parameter q 
is in fact the order parameter for the CF state. Close to 
the CF-F transition this parameter is small and one can 
expand the energy £Im/Si Eq. (^8|), in q 2 . As concerns 
the value K at the minimum, it can be found near the 
transition minimizing £Im/s at <7 = 0. As a result, the 
first terms of the expansion of the energy fl in q 2 near 
the CF-F transition can be written as 
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n^n s (K ) + n M/s (K a , q = o)- 

' F V2 ^\t\ + V2F 2 M1-K^)\t\ 



(19) 



^5/2,2 (1 - Kl 



2A 



One can check that the term proportional to q 4 is pos- 
itive, which means that the CF-F transition is of the sec- 
ond order. This is in contrast to the conclusion of Ref. 
[pf . The parameter A in Eq. (PI is 



A 



Jd 7C(3) 
/\/27^D 3 ~ 2tt 2 



(20) 



According to the Landau theory of phase transitions the 
transition from the F state (q = 0) to the CF state (q ^ 
0) should occur when the coefficient in the second-order 
term turns to zero. The phase diagram for the variables 
h and J, Eqs. ([[t], |2(]), is represented in Fig.l. The 
curves are plotted for different values of |t|. The function 
S7(g) has only one minimum at qo continuously going to 
zero as the system approaches the transition point. This 
demonstrates that the transition is of second order. Not 
close to the transition point Q ~ 

The stiffness J for materials like Fe and Ni is 
w 607^/A. Using the data for Nb T c = lOA, v F = 
1, 37.10 8 cm/s, setting I = 100A, d = lOA, and h = 10 4 K, 
which is proper for iron, and assuming that the Fermi 
velocities and energies of the ferromagnet and supercon- 
ductor are close to each other we obtain a w 25 and 
A ~ 6.1CU 3 . It is clear from Fig.l that the CF state is 
hardly possible in the Fe /Nb structure studied in |l(| . 

How can one explain the decay of the average mag- 
netic moment below T c observed in that work? This can 
be understood if one assumes that there were "islands" in 
the magnetic layers with smaller values of J and/or h. A 
reduction of these parameters in the multilayers Fe/Nb 
is not unrealistic because proximity to Nb leads to for- 
mation of non-magnetic "dead" layers j| , and can affect 
the parameters of the ferromagnetic layers, too. If the 
CF state were realized only on the islands, the average 
magnetic moment would be reduced but remain finite, 
which would correlate with the experiment [jlo|l . One can 
also imagine islands very weakly connected to the rest 
of the layer, which would lead to smaller energies of a 
non-homogeneous state. 

Another possibility to observe the CF state would be to 
use multilayers with a weaker ferromagnet. A good can- 
didate for this purpose might be Gd/Nb. The exchange 
energy h in Gd is h ps 10 3 Af and the Curie temperature 
and, hence, the stiffness J is 3 times smaller than in Fe. 
So, one can expect a « 2.5 and A « 2.10~ 3 . Using Fig.l 
we see that the CF phase is possible for these parameters. 
One can also considerably reduce the exchange energy h 
in V\~ x Fe x /V multilayers || varying the alloy composi- 
tion. Hopefully, the measurements that would allow to 
check the existence of the CF phase in these multilayers 
will be performed in the nearest future. 




a 

FIG. 1. Phase diagrams (A, a). The area above (below) the 
curves corresponds to the F (CF) state 

In conclusion, we studied a possibility of the CF state 
in (S/F) multilayers. We derived a phase diagram that 
allows to make definite predictions for real materials. 
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